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Abstract The idea that chaotic set of quantum vortices can mimic classical tur- 
bulence, or at least reproduce many main features, is currently actively being 
developed. Appreciating significance of the challenging problem of the classical 
turbulence it can be expressed that the idea to study it in terms of quantized line 

flj ' is indeed very important and may be regarded as a breakthrough. For this reason, 

r^ ^ this theory should be carefully scrutinized. One of the basic arguments supporting 

this point of view is the fact that vortex tangle decays at zero temperature, when 
the apparent mechanism of dissipation (mutual friction) is absent. Since the all 

jrt ■ possible mechanisms of dissipation of the vortex energy, discussed in the literature, 

are related to the small scales, it is natural to suggest that the Kolmogorov cascade 
takes the place with the flow of the energy, just as in the classical turbulence. In 
'■^ ' a series of recent experiment attenuation of vortex line density was observed and 

C I authors attribute this decay to the properties of the Kolmogorov turbulence. In the 

present work we discuss alternative possibility of decay of the vortex tangle, which 
is not related to dissipation at small scales. This mechanism is just the diffusive 
like spreading of the vortex tangle. We discuss a number of the key experiments, 
considering them both from the point of view of alternative explanation and of 
the theory of Kolmogorov turbulence in quantum fluids. 
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1 Introduction 



The diffusion-like decay of the vortex tangle is closely related to the hypothet- 
ical connection between the classical (Kolmogorov) and superfluid turbulence. 
The idea that chaotic set of quantum vortices can mimic classical turbulence, 
or at least reproduce many main features, is currently actively being developed 
[T],[2],[3],[l],[5],[6]. In principle, this idea had been discussed early (see e.g. famous 
textbook by Frisch [7]), as an alternative version of the problem of classical tur- 
bulence. But only now, when the new powerful experimental methods in quantum 
fluids appeared, this idea can be checked experimentally and it seems to be very 
attractive. Appreciating significance of the challenging problem of the classical 
turbulence it can be expressed that the idea to study it in terms of quantized line 
is indeed very important and may be regarded as a breakthrough. For this reason, 
this theory should be carefully scrutinized. 

One of the basic arguments supporting the idea of Kolmogorov turbulence in 
quantum fluids is the fact that vortex tangle decays at zero temperature, when 
the apparent mechanism of dissipation (mutual friction) is absent. Numerical and 
experimental observations of decay of the tangle at zero temperature are presented 
in papers [8] , [9] , [10] , [11] . The physical mechanisms of the dissipation can be var- 
ious, some approaches and ideas such as a cascade-like break-down of the loops, 
Kelvin waves cascade, acoustic radiation, reconnection loss, etc., have been dis- 
cussed in detail in recent review [12] . It is remarkable that all of these mechanisms 
are realized only on a very small scale. Therefore, it is natural to suggest that 
the Kolmogorov cascade occurs with the flow of energy, just as in the classical 
turbulence. The mentioned experimental works on decay of the vortex tangle are 
interpreted from point view of the decay of classical turbulence. In works by Skr- 
bek [2] it was developed approach, which relates the attenuation of the energy in 
classical turbulence to the decay of the vortex line density and predicts the tem- 
poral dependence of this attenuation. 

In the present work we discuss the alternative mechanism of decay of the vortex 
tangle, which is not related to dissipation at the small scales. This mechanism is 
just the diffusive spreading of the vortex tangle. We applied the diffusion equation 
obtained early ([16]) to describe the decay of superfluid turbulence in the listed 
above experiment. Our calculations enable us to conclude that mechanism of dif- 
fusion can describe correctly the attenuation of the vortex line density. Besides, 
we discuss these experiments from position of theory Kolmogorov turbulence. 
In parallel, to check our hypothesis, we perform direct numerical simulation of the 
evolution of the vortex tangle, originally concentrated in the restricted domain 
(See Fig.l). We were carefully monitoring all mechanisms of the decrease of the 
total vortex length (reconnection procedure, inserting or removing of intermediate 
points, elimination of small loops etc.). It is found that the most effective mecha- 
nism is related with the spreading of the vortex tangle, and the rate of change of 
the total length is occurring in accordance with the diffusion equation, described 
above. 
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Fig. 1 Diffusion of a vortex tangle at different times. It is clearly seen as the vortex loops 
escape from the volume carrying away the line length and energy. 



2 Diffusion Equation 

In this section we very briefly describe main ideas leading to the diffusion theory 
of the vortex loops, details can be found in paper [16]. Vortex loops composing the 
vortex tangle can move as a whole with some drift velocity V; depending on their 
structure and their length. The flux of the line length, energy, momentum etc., 
executed by the moving vortex loops takes place. In the case of inhomogeneous 
vortex tangle the net flux J of the vortex length due to the gradient of concentra- 
tion of the vortex line density C{x,t) appears. The situation here is exactly the 
same as in classical kinetic theory with the difference being that the "carriers" 
are not the point particles but the extended objects (vortex loops), which possess 
an inflnite number of degrees of freedom with very involved dynamics. We offer 
to fulfill investigation basing on the supposition that vortex loops have the Brow- 
nian or random walking structure with the generalized Wiener distribution (see 

mmm)- 

To develop the theory of the transport processes fulfilled by vortex loops (in spirit 
of classical kinetic theory) we need to calculate the drift velocity VJ and the free 
path \{l) for the loop of size I. Referring to the paper [16] we write down here the 
following result. The drift velocity Vi for the loop of size Z.is 



Vi = Cyl3/, 



(1) 



Quantity /? is {k/4:'!t)\ii{C~ ' /ag), where k is the quantum of circulation and 
ao is the core radius, C„ is numerical factor of the order of unity. The ,^o is the 
parameter of the generalized Wiener distribution, it is of order of the interline 
space £~^ ". The free path X{1) for loop of length I is: 



A(/) = l/2lbmC. 



(2) 
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Here, fom is the numerical factor , approximately equal to bm » 0.2 It is seen 
that free path X(l) is very small, it implies only very small loops give a significant 
contribution to transport processes. Knowing the averaged velocity Vj of loops, and 
the probability P{x) (both quantities are I -dependent), we can evaluate the spatial 
flux of the vortex line density £,, executed by the loops. The procedure is very 
close to the one in the classical kinetic theory, with the difference being that the 
carriers have different sizes, requiring additional integration over the loop lengths. 
Referring again to paper [16 _ we write the flux J of vortex line is proportional to 
V£ and, correspondingly, the spatial-temporal evolution of quantity C obeys the 
diffusion type equation 

^ = D.V'C, (3) 

Our approach is a fairly crude to claim a good quantitative description. However, 
if we are to adopt the data grounded on the exact solution to the Boltzmann type 
"kinetic" equation for vortex loops distribution ([15]), we conclude that Cd ~ 2.2. 
Further we use the ([3]) to describe the decay of superfiuid turbulence in various 
experiments including numerical simulations. 



3 Discussion of the experimental data and numerical simulations 

In this section we discuss several experiments and numerical simulation on the 
decay of the vortex, which are usually considered as the ground for the Kol- 
mogorov decay of the superfiuid turbulence. They are the so called Lancaster 
experiment [10|, the Manchester experiment [H] and the Osaka numerical simula- 
tion [g. 



3.1 Osaka numerical simulation 

Results of numerical simulation on the dynamics of the vortex tangle at zero tem- 
perature made in Osaka [8_ , are frequently considered as a base for conclusion that 
for zero temperature decay of the vortex tangle occurs via various mechanisms re- 
alizing at small scales. Attenuation of the vortex line density obtained in numerical 
simulation ,8J is depicted in the upper picture of the Fig. 2. 

Let us briefly analyze the situation presented in paper [8] , and demonstrate that 
none of the currently discussed mechanisms of the "homogeneous" decay of the 
vortex tangle at zero temperature, can be applied to this work. Thus the Kelvin 
waves cascade could not be a reason for the vortex tangle decay in numerical sim- 
ulation [8], simply because the space resolution A^ = 2 * 10~^ cm was too large 
to monitor the region of large wave numbers, required for observation of a gener- 
ation of higher harmonics. Similarly, the acoustic radiation could not be a reason 
for the vortex tangle decay, because the compressibility had not been included in 
the numerical scheme. As for the loss of the line length during reconnection, the 
real effect of the length loss can be obtained only on the basis of more rigorous 
theory, e.g. with the use of the Gross-Pitaevskii equation. It is known, however, 
that an artificial loss of length is possible, due to realization of the reconnection 
procedure. This effect, however, should depend on the space resolution, whereas it 
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Fig. 2 Comparison with the numerical simulation by Tsubota et.al.|8] In the upper picture 
it is depicted the attenuation of vortex line density obtained in numerical experiment |8]. In 
the lower picture it is depicted the same quantity calculated with the use of equation ((Sj 
without the auxiliary term, the diffusion constant was equal to Cd ~ 2.2 * 10~^ cm? j s. We 
calculated the two-dimensional evolution of the vortex line density in the 1 cm square resolving 
numerically equation ((3j with the boundary conditions accounting the back radiation from the 
solid walls (see for details I16| ). 



was proven that the rate of decay did not depend on it. Feynman cascade of con- 
sequent breaking down of vortex loops was imitated in [8 with elimination of very 
small loops (with sizes of few Zi^) . But the tha total amount of these events was 
is too small to describe the decay of vortex line density, observed in the numerical 
experiment. 

To resume, one concludes that none of the discussed mechanisms could be the 
reason of the "homogeneous decay" of the vortex tangle in numerical simulation[8]. 
Thus, the nature of attenuation of the vortex line density in [81 remained unclear. 
The only mechanism capable of explain the decrease of the vortex line density /I(i) 
is just spreading of the tangle and escaping of the vortex loops from the bulk of 
helium. To check our supposition, we calculated the evolution of the vortex line 
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density under conditions taking in work [8] resolving equation ^. The problem 
had been resolved numerically, the result is depicted in the lower picture of Fig. 2. 
Comparison of the upper and lower pictures of Fig. 2, enabled us to conclude that 
the diffusion spreading describes satisfactorily the evolution of the vortex tangle, 
without any additional mechanism. 

Resuming this subsection we can state that (i) there is no convincing evidence 
in favor of existence of cascade-like transfer of the vortex length (and energy)in 
direction of small scales and (ii) decay of the superfluid turbulence is quite well 
described by a diffusion mechanism. 



3.2 Lancaster experiment 

Let us now discuss the recent experiment on decay of the vortex tangle at very 
low temperatures [10] . The authors of this work reported the attenuation of vortex 
line density in superfluid turbulent helium, ^He-B. In the upper picture of Fig. 3, 
we displayed Fig 2 of work [TU] , showing results of measurements on the temporal 
behavior of the average vortex line density C{t), (solid curves, see [1^ for details). 
Authors compare their data with the theory by Skrbek |2J and conclude that decay 
of the vortex tangle occurs in accordance with theory of classical (Kolmogorov) 
turbulence. The main argument is that the long time attenuation of the vortex 
length is described by the same line A ("limiting behavior") with the t~^'^ depen- 
dence (almost independently on initial level). The according kinematic viscosity is 
v' ~ 0.2k. 

It is necessary to note the following circumstance. The one of the standard vision 
how the set of the vortex filament can imitate the classical turbulence is that the 
lines are unified into the bundles (containing many lines) . The set of many bun- 
dles of various sizes randomly move, imitating the dynamics of eddies in classical 
turbulence. Other view is that in the dense vortex tangle there is polarization 
" indistinguishable by glance" and these polarized vortices also reproduce the eddy 
dynamics. It is necessary to note the following circumstance. The one of the stan- 
dard vision how the set of the vortex filament can imitate the classical turbulence 
is that the lines are unified into the bundles (containing many lines). The set of 
many bundles of various sizes randomly move, imitating the dynamics of eddies 
in the classical turbulence. Other view is that in the dense vortex tangle there is 
an averaged partial polarization of lines "indistinguishable by glance" and these 
polarized vortices also reproduce the eddy dynamics. Let us consider the situation 
in the Lancaster experiment more attentively. Let us take some "intermediate" 
value of the vortex line density C{t) = 0.5* 10^ cm~^, where all lines are collapsed 
into the "limiting" , universal behavior (line "A" in the upper picture of Fig. 3). For 
this value the interline space is equal to C~ " « 4.5 x 10~^ cm. The latter implies 
that in the volume with size d = 1.5 mm (offerred by authors as the region, where 
the vortex tangle evolves) we have approximately d/L~^''^ w 3 lines. Of course, 
this amount is not enough to form many bundles with very "dense array of vortex 
lines" , which are necessary to "mimic classical turbulence" . On the same ground 
one can assert that it is hardly possible to say about the partial polarization of 
lines in the dense tangle " indistinguishable by glance" . 

Another fact is that about one third of the "limiting line A" is occupied by the 
developed fiuctuations, which can blur the true dependence. 
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Fig. 3 Comparison with experiment 1101 . See the text Comparison with experiment 1101 . In 
the upper picture we displayed Fig 2 of work |10| . showing results of measurements on the 
temporal behavior of the average vortex line density £(t), (solid curves, see 1101 for details). 
We calculated the same quantity resolving the diffusion equation ^ , with the use the boundary 
condition, which corresponds to the smearing of the vortex tangle into ambient space. It is 
known that for this condition the solution of three-dimensional problem is just production of 
three one-dimensional solutions. The straight line in the lower picture exactly corresponds to 
line A, in the upper picture (this line was named by the authors of paper 1101 as a "limiting 
behavior"). 



Consequently supposing that the decay is reahzed by the escaping of the vortex 
loops we had estimated the contribution into attenuation of the vortex line density, 
due to the pure diffusion mechanism. We use the classical solution (see for details 
|16p. Results (with comments) are depicted in Fig. 3. We again can conclude that 
the diffusion spreading describes satisfactorily the evolution of the vortex tangle, 
without any additional mechanism. 

Concluding this subsection we again can state that (i) Interpretation of experi- 
ment is not fully consistent and cannot serve as convincing evidence in favor of 
existence of cascade-like transfer of the vortex length (and energy)in direction of 
small scales, and (2)decay of the superfluid turbulence is quite well described by 
the diffusion type mechanism. 
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Fig. 4 Comparison with experiment |11| . See the text Comparison with experiment [6] . In the 
left picture of Fig. 6, it is depicted the temporal behavior of the average vortex line density 
Cav (t) is depicted. We calculated the same dependence on the basis of the diffusion equation 
SSl . with the boundary boundary accounting back radiation of loops from the solid walls. The 
full three dimensional problem had been resolved numerically. 



3.3 Manchester experiment 



The next experiment which we would hke to discuss and which is also frequently 
referred as the evidence of the Kolmogorov turbulence imitated by quanitized vor- 
tex lines is the work [IL. In this work the decay of vortex tangle in He-II was 
observed in the closed cube with soUd walls. Results are collected in the upper 
picture of Fig.4,where the temporal behavior of the average vortex line density 
Cav{t) is depicted. Authors compare their data with the theory by Skrbek [2] and 
conclude that decay of the vortex tangle occurs in accordance with theory of clas- 
sical (Kolmogorov) turbulence. The according kinematic viscosity is i/ '^ 0.003k , 
which about two orders smaller, than obtained in the Lancaster experiment. Au- 
thors of [11] assumed that the source of this discrepancy is that the turbulence 
observed in [10^ is not homogeneous, and the size of the energy-containing eddies 
may differ from the spatial extent of the turbulence, so that the value of v' ob- 
tained in '10] was ambiguous. 

Coming to the previous subsection we can say that situation with the number of 
vortices is better than in the Lancaster experiment. It is should be noted however 
that it is merely due to larger size of the experimental cell. If to take parts of 
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experimental cell with sizes of the order 3 mm we again meet the of very dilute 
vortex tangle. Another remark is related to the Volovik observation ^7\. At the 
low temperature in '^He-B and almost at any temperature in *He, there should the 
so called Vinen (not Kolmogorov) turbulence, and, the whole ideology grounded 
on the Skrbek theory is not applicable in the case of Manchester experiment. 

In our opinion the large differnce in the value of kinematic viscosity is v' in 
works [TT] and \Wl (In fact it is just difference in the total time of the decay) 
is related to (i) the different sizes of the cells, where the superfluid turbulence is 
activated, and (ii) to the different boundary conditions (solid walls in [11], and 
absence of boundaries in [TU]). These two facts definitely point out that the decay 
has the diffusive nature. We calculated the decay of the vortex tangle on the basis 
of the diffusion equation ([3]). The fully three-dimensional problem was resolved 
numerically (see for details [16|). the result shown in the lower picture of Fig. 4. 
It can be seen that the decay of the vortex tangle, due to diffusion, describes both 
quantitatively and qualitatively the features observed in the experiment. First of 
all, the whole qualitative behavior of lines agrees with diffusive-like attenuation. 
In particular, there is a plateau, which is changed with the fast decay of the 
tangle. Full decay of the superfiuid turbulence occurs in times, which are in a very 
good agreement, predicted on the basis of the diffusion approach elaborated here. 
The slope of the curve in the interval of the most intensive decrease, shows the 
dependence close to ~ i~'^", which is also typical for diffusion. 
Resuming this subsection we again claim that (i) experiment is not fully self- 
consistent and (ii) the diffusion mechanism describes well the attenuation of the 
vortex line density. 



3.4 Novosibirsk numerical simulation 

We fulfilled the direct numerical simulations of the evolution of the nonuniform 
vortex tangle, originally concentrated in the restricted domain, at zero tempera- 
ture. The numerical simulation is performed with the use of local induction approx- 
imation (LIA). An algorithm, which is based on consideration of crossing lines, is 
used for vortex reconnection processes. The calculations are carried out in an infi- 
nite volume. Result of the numerical simulations is depicted in Fig.l, it is clearly 
seen how the vortex loops escape from the volume and vortex line density dimin- 
ishes inside the bulk. We properly studied the change of the total vortex length 
in the initial domain. In particular we were carefully monitoring all mechanisms 
leading to the loss of the length, such as the change of length owing to reconnec- 
tion processes, the eliminations of very small loos below the space resolution, the 
change of length due insertion and removing of points to supply numerical algo- 
rithm stability and etc. Results of this monitoring are depicted in Fig. 5. From this 
figure we conclude, that the master mechanism of decrease of the vortex length 
inside initial volume is related to the escaping of vortex loops, with is realized in 
the diffusion like manner, and the other mechanisms listed above do not affect 
appreciably on the spreading of the vortex tangle. 

We compared results of our numerical simulations with the theory of diffusion of 
the vortex tangle described above. We obtained, that the evolution of of the length 
in initial domain is satisfactory described by the with the diffusion equation ^. 
Result of comparison is depicted in Fig 6. The good agreement with the exper- 
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Fig. 5 Contribution of various mechanisms into decrease of total length. 1-Dccrease of total 
length, 2-Decrease of total length inside domain (spherical domain with the radius Rj 0.008 cm) 
due to escaping of the small loops, 3-Decrease of total length due to elimination of small loops 
(3 mesh sizes ), 4-Change of total length due to artificial procedure - inserting or removing 
of intermediate points, 5-Change of total length due to reconnection procedure, 6 -Change of 
total length due to motion after reconnection, 7- Total length inside domain 
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Fig. 6 Total length l{t) = j C{r,t)dr inside domain obtained in numerical simulation (left 
picture) and calculated on the base of the diffusion equation (for spherical domain) 



imental data and numerical simulation enables us to conclude that the diffusion 
process plays a dominant role in the free decay of the vortex tangle in the absence 
of the normal component. 



4 Conclusion 



Our main conclusion would be formulated as follows. The both experimental and 
numerical data on decay of the superfluid turbulence discussed in the literature 
cannot be regarded as the firm evidence for the classical turbulence mechanism, ac- 
companied by the Kolmogorov cascade of the energy to region of very small scales. 
As we demonstrated, they are not fully self-consistent even in the frame of the ac- 
cepted approach. At the time we had shown that all experimental data on decay of 
the vortex tangle agree well with the diffusion mechanism without any additional 
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assumptions. Our numerical results confirm this point of view, demonstrating that 
all possible mechanisms of the loss of the vortex length (and, correspondingly, the 
vortex energy) are considerably less than the loss of energy due to escaping of the 
loops from the volume. 

Acknowledgements Authors are grateful to participants of workshop " SYMPOSIA ON SU- 
PERFLUIDS UNDER ROTATION" for the very useful discussion. 
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